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ABSTRACT 


This  report  is  a review  of  evidence  for  anomalous  attenuation  of  short- 
period  seismic  body-waves  In  the  upper  mantle  under  the  Western  United  States. 
The  attenuation  reduces  the  amplitudes  of  these  waves  by  a factor  of  about  two 
relative  to  those  observed  In  stable,  continental  shield  type  areas.  The  mag- 
nitude of  this  amplitude  anomaly  Is  consistent  with  spectral  differences  of 
body-waves  observed  In  the  same  region,  assuming  a simple  exponential  atten- 
uation law  with  a constant  quality  factor  Q.  The  most  likely  cause  of  the 
attenuation  Is  partial  melting  In  the  upper  mantle  under  the  Western  United 
States.  This  hypothesis  Is  supported  by  measurements  of  travel-time  delays, 
upper  mantle  conductivity  and  heat  flow.  The  report  also  Includes  a consider- 
able number  of  new  results  obtained  recently  at  the  Seismic  Data  Analysis 
Center. 


1 


TABLE  OF  CONTENTS 


ABSTRACT 
LIST  OF  FIGURES 
LIST  OF  TABLES 
INTRODUCTION 

DIRECT  EVIDENCE  FOR  THE  EXISTENCE  OF  MAGNITUDE  BIAS  FOR 
EVENTS  IN  THE  WESTERN  UNITED  STATES  (WUS)  VS.  EVENTS  IN 
CONTINENTAL  SHIELD  AREAS  (EUS  IN  PARTICUIJUl) 

SHIFT  OF  POPULATIONS 

SUPPORTING  EVIDENCE  FOR  HIGH  ATTENUATION  UNDER  THE  WUS 
General 

Attenuation  of  Long-Period  Body  and  Surface  Waves 
Travel  Time  Delays 
Velocity  Structure 

Upper  Mantle  Resistivity  and  Heat  Flow 
GENERAL  COMMENTS 
ACKNOWLEDGEMENTS 


Page 

2 

4 

6 

7 

9 

36 

38 

38 

38 

42 

46 

53 

56 

58 


REFERENCES 


59 


LIST  OF  FIGURES 


Figure  Title 

1 Short-period  magnitude  residuals  in  the  United  States 
based  on  LRSM  bulletins  and  compiled  by  Booth,  liarshall 
and  Young,  1974. 

2 Summary  of  9“  < A < 35*  t*  values  classified  with  re- 
spect to  EUS  and  mixed  path  types.  Der  and  McElfresh 
(1976). 

3 Reduced  magnitudes  from  explosions  GNOME  and  SALMON  as 
functions  of  t*.  Solid  line  corresponds  to  the  ampli- 
tude factor  exp(-ii)t*/2)  for  f = 1 Hz.  Der  and  McElfresh 
(1976b).  Dashed  line  gives  the  theoretically  calculated 
log(A/T)  with  A corrected  for  LRSM  system  response  at 
period  T for  a 5 kiloton  explosion  in  salt  using  the  re- 
duced displacement  potential  of  Gnome  from  Werth  and 
Herbst  (1963)  and  the  calculation  techniques  from  Bland- 
ford  (1976). 

4 Summary  of  t*  values  for  A > 35°  listed  in  tables  II 
and  111. 


Relative  plotted  against  relative  t*  for  five  Novaya 
Zemlya  events. 

LTM  discriminants  at  NORSAR  for  events  located  in  various 
regions  of  the  world  (after  Noponen  1975).  Explosion 
and  earthquake  observations  are  plotted  with  filled  and 
open  symbols,  respectively.  (a)  Events  in  southwestern 
North  America.  (b)  Explosions  in  western  Russia,  earth- 
quakes in  southern  Greece  (circles)  and  Turkey  (tri- 
angles). (c)  Events  in  Central  Asia.  (d)  Earthquakes 
at  Kuriles  (circles)  and  Sakhalin  (triangles).  Almost 
all  of  the  events  in  the  figure  are  shallow,  and  none  are 
deeper  than  100  km.  Lines  denote  LTM  theoretical  values 


computed  for  t*  values  indicated  in  the  figure  for  explo- 
sions (solid  lines)  and  earthquakes  (dash  lines). 


6b  Summary  of  all  t*  data  in  Table  I as  functions  of  epi-  30 

central  distance. 

7 Map  of  short-period  S-wave  attenuation  after  Der,  Mass6  32 

and  Gurskl  (1975). 

8a  The  distribution  of  S-wave  periods  in  the  WUS  and  EUS.  33 

(Average  of  periods  for  each  event  is  subtracted) . 

8b  Two  Q models  for  short-period  P-waves  which  satisfy  the  35 

data  presented  in  this  report. 


-4- 


LIST  OF  FIGURES  (Continued) 


I 

I 

* 

' i 


1 

1 • A 

1 1 

Figure  No 

. Title 

Page 

V 

9 

relationships  for  earthquakes  and  explosions  in 
various  regions  of  the  world  (after  Basham  1969), 

37 

h 

♦ \ 

V 

10 

t*  values  for  long-period  P-waves  in  the  United  States 
(after  Solomon  and  ToksHz  1970). 

39 

11 

t*  values  for  long-period  S-waves  in  the  United  States 
(after  Solomon  and  ToksHz  1970). 

39 

■‘i 

12 

Comparison  of  S and  P differential  attenuation  at  U.S. 
stations.  Vertical  and  horizontal  error  bars  delimit 
the  range  of  calculated  values  for  each  station.  Cor- 
relation coefficient  is  0.67  (after  Solomon  and  ToksHz 
1970). 

40 

♦ 

13 

Limits  of  100/Q  values  as  functions  of  depth  for  three 
major  regions  or  the  United  States  (after  Lee  and  Solo- 
mon 1975). 

43 

1 

1 

14 

Travel-time  residuals  for  short-period  P-waves  from 
deep  earthquakes  at  various  stations  in  the  United 

States  (after  Julian  and  Sengupta  1973). 

44 

1 

15 

S-wave  travel-time  residuals  at  various  stations  in  the 
United  States  (after  Doyle  and  Hales  1967), 

44 

16 

S and  SKS  travel-time  residuals  in  the  United  States 
(after  Hales  and  Roberts  1970). 

45 

17 

P and  S-wave  travel-time  delays  indicating  a linear  cor- 
relation between  the  two  (Doyle  and  Hales  1976). 

47 

18 

Selected  P-wave  velocity  models  for  the  Western  United 
States. 

48 

19 

Selected  P-wave  velocity  models  for  the  Eastern  United 
States. 

50 

- 

20 

P velocities  under  the  United  States, 
n 

52 

♦ 

21 

Upper  mantle  conductivity  map.  Density  of  stippling 
proportional  to  conductivity. 

54 

22 

Heat  flow  measurements  in  the  Western  United  States. 

55 

A 


LIST  OF  TABLES 

Table  No.  Page 

I Summary  of  t*  measurements.  12 

Ia,b  t*  values  for  9 < A < 35°  displayed  in  Figure  2.  15 

II  t*  values  for  shield  and  stable  platform  paths  23 

A > 35°. 

Ill  t*  values  for  mixed  (shield  and  stable  platform  to  24 

Basin  and  Range)  paths  A > 35°. 

IV  t*  values  for  paths  to  the  Colorado  Plateau.  25 


-6- 


INTRODUCTION 

Many  factors  influence  amplitudes  of  seismic  waves  from  explosions  and 
earthquakes.  The  identification  of  an  explosion  and  estimation  of  its  yield 
involve  amplitude  measurements  on  at  least  two  kinds  of  waves  (the  P-wave  ar- 
rival and  the  Rayleigh  wave  amplitude) . In  view  of  the  recently  proposed 
Threshold  Test  Ban  Treaty  (TTBT) , which  limits  the  yields  of  nuclear  explo- 
sions, it  is  important  to  know  how  amplitudes  of  seismic  waves  depend  on  the 
region  in  which  they  were  fired.  There  is  a considerable  amount  of  evidence 
that  indicates  that  the  amplitudes  of  waves  at  the  teleseismic  distance  de- 
pend strongly  on  the  geophysical  properties  of  the  source  region.  To  esti- 
mate reliably  the  yield  of  explosions  one  must  know  the  magnitude  of  such 
anomalies  for  each  source  region. 

The  body  wave  magnitude  (m^)  is  derived  by  the  well-knoxm  formula 

m^  = log^Q(A./T)  + B(A) 

where  A is  the  amplitude  of  the  P-wave,  T is  its  dominant  period  and  B is  a 
distance  dependent,  empirical  term  which  approximately  corrects  for  geometri- 
cal spreading  and  attenuation  in  an  average  Earth.  All  the  quantities  used 
in  this  formula  strongly  depend  on  transmission  factors,  many  of  which  cannot 
be  easily  estimated,  such  as  the  velocity  structure,  variations  in  attenua- 
tion, receiver  structure,  etc.  As  a consequence,  raw  m^  values  for  any  given 
event  show  a large  scatter. 

Not  all  of  the  factors  influencing  m^^  for  a given  event  will  be  dis- 
cussed in  this  report.  The  effects  of  near  source  medium,  device  yield, 
source  asymmetry,  surface  reflection  or  spall  and  near-surface  geology  are 
excluded  from  this  discussion. 

We  also  exclude  the  effects  of  multipathing  and  scattering  from  this 
discussion  on  the  basis  that  such  effects  are  essentially  uncontrollable  for 
single  sites  and  extremely  difficult  to  evaluate  for  arrays.  One  can  assume 
that  the  scatter  due  to  such  effects  can  be  reduced  by  taking  averages  of  m^ 
over  many  stations.  In  this  report  we  will  examine  broad  regional  differ- 
ences in  m^  values.  Special  attention  will  be  given  to  the  properties  of  the 
mantle  under  Nevada  Test  Site  (NTS)  as  compared  to  other  nuclear  test  sites. 
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An  extensive  review  of  the  magnitude  problem  was  given  recently  by  Rodean 
(1976)  who  summarized  all  the  relevant  research  to  date.  This  report  will 
show  much  of  the  Illustrative  material  from  the  sources  quoted  by  Rodean,  and 
will  discuss  some  relationships  among  various  kinds  of  data  as  well  as  some 
discrepancies.  Beyond  qualitative  geographical  correlations  between  the 
various  relevant  geophysical  parameters  It  will  be  shown  that  the  observed 
P-wave  amplitude  variations  are  consistent  with  variations  of  Q In  the  upper 
mantle  as  deduced  from  spectral  measurements  of  P-waves. 


Rodean,  H.  C.,  1976,  Seismic  yield  verification  and  a regional  vs.  m^ 
anomaly,  Lawrence  Livermore  Laboratory,  UCID-17006. 


DIRECT  EVIDENCE  FOR  THE  EXISTENCE  OF  MAGNITUDE  VARIATION  FOR  EVENTS 
IN  WESTERN  UNITED  STATES  (WUS)  VS.  EVENTS  IN  CONTINENTAL  SHIELD  AREAS 

(EUS  IN  PARTICULAR) 


Many  investigators  have  observed  that  the  body-wave  magnitudes  for  tele- 
k seismic  events  measured  in  the  Western  United  States  (WUS)  differ  from  those 

measured  in  the  Eastern  United  States  (EUS)  (Evernden  and  Clark,  1970;  Booth, 
Marshall  and  Young,  1974;  North,  1976).  Body-wave  magnitude  measurements  in 
the  WUS  were  consistently  0.25  magnitude  units  lower  than  in  the  EUS.  The 
average  regional  differences  computed  by  other  authors  are  in  the  range  of 
0.25-0.4  magnitude  units.  Figure  1 shows  the  magnitude  terms  for  short-per- 
iod P-waves  at  selected  LRSM  stations  as  determined  from  about  400  events  by 
Booth,  Marshall,  and  Young  (1974).  Although  magnitude  variations  can  have 
many  causes,  such  as  thick  low-velocity  sedimentary  surface  layers  at  the 
given  stations,  or  focusing  or  defocusing  of  rays  by  local  inhomogeneities, 
we  shall  see  that  the  most  plausible  explanation  of  the  general  trend  of 
short-period  P-wave  magnitude  variations  and  associated  phenomena  in  the  WUS 
^ is  anelastlc  attenuation  due  to  a viscous,  partially  molten,  low-velocity 

* zone  (LVZ)  under  the  WUS. 

Anelastic  attenuation  decreases  the  amplitude  of  a given  frequency  com- 

* ponent  of  a wave  by  the  factor; 

exp  (-TTf  t*) 

where  f is  frequency,  t*  = t/Q  where  t is  the  travel  time  of  a seismic  wave 
and  Q is  the  average  quality  factor  along  a given  seismic  ray.  Since  Solo- 
mons  and  ToksBz  (1970)  obtained  lower  Q values  for  long  period  waves  than  those 
deduced  in  this  report  the  quality  factor  Q may  be  a function  of  frequency, 
but  we  shall  assume  it  to  be  a constai'.t  between  frequencies  .5  and  4 Hz,  the 
frequency  band  of  short-period  instruments  used  for  measuring  m^.  An  effective 

Evernden,  J.  and  D.  M.  Clark,  1970,  Study  of  teleseismic  P.  II.  Amplitude  data, 
■ Phys.  Earth.  Planet.  Int. , v.  4,  p.  24-31. 

Booth,  D.  C.,  P.  D.  Marshall,  and  J.  B.  Young,  1974,  Long  and  short  period 
I amplitudes  from  earthquakes  in  the  range  0'’-li4°,  Geophys . U.  Astr . 


Soc. , V.  39,  p.  523-538. 

North,  R.  G.  (1976)  Station  biases  in  body  wave  magnitude.  (manuscript). 
Solomon,  S.  C.,  and  M.  N.  ToksHz,  1970,  Lateral  variation  of  attenuation  of 


I 

I 

i 

i 

i 


I 


value  of  t*  can  be  measured  directly  by  comparing  the  observed  spectra  of  P- 
waves  with  that  of  the  source,  derived  from  close-in  measurements  (RDP).  Re- 
lative differences  of  effective  t*  can  be  found  by  comparing  the  spectra  at 
various  stations. 

If  the  observed  magnitude  variation  is  due  to  anelastic  attenuation  it 
has  to  be  consistent  with  the  t*  values  derived  from  spectral  ratio  measure- 
ments. A hot,  partially  molten  upper  mantle  will  lower  the  shear  modulus  of 
the  material  relative  to  a solidified  cold  mantle  and  lead  to  anelastic  energy 
losses  in  shear  deformation.  This  condition  will  reduce  the  S-wave  amplitudes 
and  velocities  even  more  drastically  than  those  of  P-waves  in  the  same  fre- 
quency range.  The  values  of  anelastic  t*  for  S-waves  should  be  three  to  four 
times  those  for  P-waves  (Solomon  and  Toks8z,  1970).  Anelastic  attenuation  as 
a cause  of  amplitude  variation  will  also  imply  reciprocity,  so  that  sources 
located  in  the  WUS  will  seem  to  a distant  observer  to  release  less  energy. 

In  Table  I we  show  all  estimates  of  effective  t*  values  known  to  us  at 
the  present  time.  We  include  values  from  the  literature  as  well  as  from  our 
own  research.  In  some  cases  we  estimate  t*  from  data  or  Q models  published 
where  t*  was  not  stated  explicitly.  We  also  indicate  the  computational  me- 
thods and  assumptions  used  to  derive  t*,  the  epicentral  distance  range,  and 
the  nature  of  paths  and  sources.  The  model  of  Archambeau  et  al.  (1969a, b) 
shows  Q values  which  are  too  high  relative  to  those  derived  by  us  in  the  Basin 
and  Range,  since  the  observed  amplitude  loss  was  interpreted  in  terms  of  a 
horizontally  homogeneous  model  where  the  same  Q model  is  assumed  at  source  and 
receiver.  The  receivers  in  Archambeau ’s  data  are  mostly  in  the  high-Q  EUS. 

In  our  interpretation  of  a causative  LVZ  practically  all  of  the  attenuation 
observed  would  be  associated  with  a value  only  half  of  that  computed  by 
Archambeau  et  al.  from  a horizontally  homgeneous  model.  The  Q along  the  as- 
cending paths  reamins  high.  (The  t*  derived  from  a complete  path  through  the 
Q model  of  Archambeau  et  al.  is  however  a correct  measure  of  the  attenuation.) 

The  t*  values  in  Table  I range  between  -0.05  and  1.  Shield  type  paths 
show  lower  values,  while  paths  including  the  Basin  and  Range  (BR)  or  tectonic 


Archambeau,  C.  B.,  E.  A.  Flinn,  and  D.  H.  Lambert,  1969a,  Fine  structure  of 
the  upper  mantle,  Geophys.  Res. , v.  74,  p.  5825-5865. 

Archambeau,  C.  B.,  E.  A.  Flinn,  and  D.  G.  Lambert,  1969b,  Detection,  analysis 
and  interpretation  of  teleseismic  signals,  1.  Compressional  waves  from 
the  SALMON  event,  Geophys.  Res. , v.  73,  p.  3877-3883. 
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NTS  Nuclear  Explosions  NTS-SDCS  Stations  BR-S  1.00  Matching  amplitude  and  Systems,  Science  and 

waveforms  of  synthetic  Software,  Inc.  Reports 

seismograms  to  observa- 
tions 
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TABLE  la,b 


WUS-EUS 

\ 

A“ 

1 

m 

NTS  - 

FKCO 

9.7 

.25 

NTS  - 

RKON 

21.2 

.49 

NTS  - 

WNSD 

13.5 

.40 

NTS  - 

KCMO 

16.5 

.59 

NTS  - 

RKON 

21.0 

.58 

NTS  - 

CPO 

24.7 

.67 

NTS  - 

FNWV 

28.9 

.66 

1 

N.  Mexico  - 

KNUT 

9.4 

.48 

- 

N.M  - 

CPCL 

10.6 

.25 

1 

N.M  - 

HLID 

14.0 

.48 

\ 

Miss.  - 

DRCO 

16.3 

.77 

Miss.  - 

IJBO 

18.6 

.59 

Miss.  - 

KNUT 

20.1 

.64 

Miss.  - 

EKNU 

22.8 

.44 

17.7 

.52 

14  values 


KNICKERBOCKER 


PILEDRIVER 


GNOME 


SALMON 


Table  Ia,b  (Continued) 


EUS 

A® 

t* 

N.M. 

- 

MFAR 

9.3 

.09 

N.M. 

- 

CUAR 

10.3 

.10 

N.M. 

- 

WKAR 

11.1 

.09 

N.M. 

- 

WNSD 

11.3 

.08 

N.M. 

- 

JSTN 

13.1 

.23 

► 

GNOME 

N.M. 

- 

CVTN 

14.1 

.16 

N.M. 

- 

NGWS 

18.1 

.25 

N.M. 

- 

BLWV 

19.3 

.24 

N.M. 

- 

DHNY 

25.0 

.28  , 

Miss. 

- 

BLWV 

9.5 

.00 

Miss. 

- 

VDIO 

11.3 

.06 

Miss. 

- 

BRPA 

12.4 

.03 

Miss. 

- 

WFMN 

12.8 

.35 

Miss. 

- 

RTNM 

13.5 

.20 

► 

SALMON 

Miss. 

- 

DHNY 

16.1 

.12 

Miss. 

- 

GPMN 

16.8 

.11 

Miss. 

- 

RKON 

19.1 

.23 

Miss. 

- 

SVQB 

28.7 

.23  , 

Avg. 

15.1 

.158 

18  values 


(T)  path  segments  are  higher.  We  purposely  excluded  from  Table  1 t*  Cor  Q) 
values  published  for  paths  Including  Island  arc  or  oceanic  ridge  segments. 
Instead,  we  center  our  discussion  on  relative  differences  between  the  WUS  and 
the  EUS  or,  in  a broader  sense,  between  pure  shield  (stable  platform)  type 
paths  and  paths  which  Include  a segment  of  the  Basin  and  Range  province. 

We  have  computed  a large  number  of  effective  t*  values  by  dividing  the 
spectra  of  the  first  nine  seconds  of  the  P-wave  signal  at  various  LRSM  sta- 
tions with  estimated  explosion  source  spectra  which  were  derived  from  scaling 
of  close-in  seismic  measurements  (Von  Seggern  and  Blandford,  1972).  These  t* 
values  compare  well  with  the  other  t*  values  published  or  implied  by  other 
authors  in  Table  I.  The  scatter  is  not  great  considering  that  we  are  dealing 
with  short-period  data,  and  that  the  t*  values  listed  were  derived  by  various 
methods. 

The  t*  values  for  the  9°  < A < 36°  distance  range  are  more  likely  to  be 

influenced  by  differences  in  the  crust-upper  mantle  velocity  structure  due  to 
travel  time  triplications  and  crustal  guided  waves  than  those  for  greater  dis- 
tances. It  is  unlikely,  however,  that  such  influences  are  very  important  for 
relative  t*  estimates.  What  we  are  measuring,  in  effect,  is  the  rate  of  fall- 
off  of  the  spectra  at  high  frequencies.  Travel  time  triplications,  crustal 
response,  multlpathlng  and  other  often  quoted  factors  may  modulate  the  spectrum 
by  introducing  nulls  and  maxima,  especially  at  shorter  epicentral  distances 
(the  values  in  Table  I with  9°  < A < 36°),  but  they  cannot  explain  the  uniform 

removal  of  high-frequency  energy  within  the  seismic  pass  band.  Seismograms  at 
the  same  epicentral  distances  have  similar  envelopes  due  to  the  same  upper  man- 
tle triplications  of  travel  times  associated  with  the  steep  gradients  at  the 
400  and  600  km  depths  that  are  identical  for  the  WUS  and  EUS  reference,  but 
the  frequency  content  of  the  seismograms  is  quite  different  due  to  the  differ- 
ent attenuations  along  the  various  paths. 

The  effect  is  visibly  quite  dramatic  on  seismograms  for  small  explosions 
where  the  source  spectrum  is  nearly  flat  within  the  frequency  range  studied 
(the  SALMON  and  GNOME  nuclear  explosions,  for  example).  For  these  events  the 
3-4  Hz  energy  is  almost  completely  absent  at  stations  with  A > 9°  along  EUS- 
WUS  paths,  while  it  is  conspicuously  present  along  pure  EUS  (shield)  type 
paths.  (Der  and  McElfresh,  1976a)  The  effect  is  still  visible  to  the  eye  but 
less  obvious  for  larger  explosions  where  most  of  the  energy  is  at  frequencies 
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f < 1.5  Hz,  but  spectral  analysis  shows  that  the  falloff  rate  at  higher  fre- 
quencies Is  considerably  steeper  for  EUS-WUS  type  paths  In  all  cases  (Der  and 
McElfresh,  1976b).  , 

Because  of  the  method  of  analysis  employed  here,  the  effects  of  surface 
reflections  are  also  small.  For  each  explosion  used  In  determining  t*  we 
assessed  the  possible  effect  of  a linear  superposition  of  pP;  In  all  the  cases 
we  analyzed.  It  seemed  to  be  minor  for  at  least  one  of  the  following  reasons: 

a)  There  are  no  obvious  nulls  in  the  spectra  which  can  be  associated 
with  pP  (most  common  case),  or  the  surface  reflection  coefficient  is 
small; 

b)  The  spectral  ratio  straight-line  fit  covers  several  nulls  and  maxima, 
thus  smoothing  out  the  pP  effect; 

c)  The  fit  is  over  one  maximum  (which  should  give  zero  slope  with  small 
attenuation  and  a flat  source  spectrum),  but  lobes  at  higher  frequen- 
cies are  missing,  indicating  high  t*;  and 

d)  We  determined  relative  t*  values  in  which  case  the  source  spectrum 
and  its  modulation  are  irrelevant. 

Next  we  will  show  that  the  observed  magnitude  variations  are  consistent  ^ 

with  the  differences  in  effective  t*  as  determined  from  spectral  ratios.  We 
first  discuss  t*  values  determined  by  SDAC  for  sources  between  the  eplcentral  * 

distance  9“  < A < 35“.  We  classify  all  t*  values  from  North  America  into  two 
groups:  "EUS"  type  paths  which  are  entirely  under  the  eastern,  tectonically 
stable  part  of  the  continent,  and  "mixed"  EUS-WUS  paths  which  cross  into  the 
WUS.  Histograms  constructed  for  these  two  groups  are  shown  in  Figure  2. 
is  a clear  separation  between  the  t*  values  of  the  two  groups.  The  mean  t* 
for  EUS-type  paths  is  0.16,  and  for  mixed  paths  it  is  0.52.  The  data  in  both 
groups  cover  the  same  distance-travel  time  range. 

Figure  3 shows  a plot  of  relative  magnitudes  (observed  m^  - log (yield)) 
of  GNOME  and  SALMON  at  various  LRSM  stations  vs.  t*  at  these  stations.  Magni- 


tude is  proportional  to  log(A/T).  (Blandford,  1976;  and  Werth  and  Herbst, 
1963)  There  is  a considerable  amount  of  scatter,  but  the  amplitudes  tend  to 


Blandford,  R.  R.,  1976,  Experimental  determination  of  scaling  laws  for  con- 
tained and  cratering  explosions,  SDAC-TR-3,  Teledyne  Geotech,  Alexan- 
dria, Virginia.  ADA  03635 


w 


log  A+B-log  (yieM) 


Figure  3. 


Reduced  magnitudes  from  explosions  GNOME  and  SALMON  as  functions 
of  t*.  Solid  line  corresponds  to  the  amplitude  factor  exp(-a)t*/2) 
for  f = 1 Hz.  Der  and  McElfresh  (1976h).  Dashed  line  gives  the 
theoretically  calculated  log(A/T)  with  A corrected  for  LRSM  system 
response  at  period  T for  a 5 kiloton  explosion  in  salt  using  the 
reduced  displacement  potential  of  GNOME  from  Werth  and  Herbst  (1963) 
and  the  calculation  techniques  from  Blandford  (1976). 


decrease  with  increasing  t*.  The  slope  of  the  solid  line  drawn  corresponds  to 
exp(-iit*),  the  attenuation  factor  at  f = 1,  the  actual  falloff  with  t*  is 
greater  because  the  dominant  frequency  of  P-wave  is  greater  than  1 Hz  for 
these  events.  The  dashed  line  shows  this  in  more  detail.  It  is  calculated  by 
the  methods  of  Blandford  (1976).  Theoretical  waveforms  are  calculated  and 
the  maximum  peak-to-peak  amplitude  is  measured.  The  interval  between  zero- 
crossings  on  either  side  of  the  maximum  peak-to-peak  excursion  is  taken  as 
the  period  T,  and  A is  then  the  measured  amplitude  corrected  for  the  system 
response  at  period  T.  (Both  lines  have  been  shifted  vertically  to  fit  the 
observed  absolute  level  of  the  data  visually.)  This  was  done  because  we 
have  not  performed  the  complete  calculations  necessary  to  predict  the  abso- 
lute level  of  the  signals  from  a specified  yield.  Regression  analysis  shows 
that  slope  zero  can  be  rejected  with  > 99%  confidence  and  that  about  half  of 
the  observed  variance  in  m^  is  due  to  t*. 

The  outlying  points  on  Figure  3 do  not  necessarily  represent  large  fluc- 
tuations in  m alone.  Errors  in  t*  measurement  also  contribute  to  the  scatter, 
b 

Figure  A shows  histograms  of  t*  values  at  teleseismic  distance  A > 35°, 

The  t*  values  used  in  this  histogram  are  summarized  in  Tables  II  and  III.  We 
have  included  all  published  t*  values  known  to  us  as  well  as  our  own.  Rela- 
tive t*  values  for  Novaya  Zemlya  events  (Novaya  Zemlya  is  known  to  be  on  a 
stable  platform  (Molnar  and  Oliver,  1969)  were  also  included  by  adding  a con- 
stant to  all  relative  values  to  make  the  average  t*  of  shield  paths  equal. 

We  subdivided  the  associated  wave  paths  into  two  groups:  the  paths  which 
start  and  end  in  tectonically  stable  regions  (shield  in  an  extended  sense) 
and  those  that  have  one  end  (ascending  or  descending)  in  the  Basin  and 
Range  province  while  the  other  end  of  the  raypath  is  in  a shield  region. 

The  difference  in  means  is  0.31.  In  the  average  of  the  shield  to  WUS  paths 
we  included  two  extremely  high  t*  values,  the  Mast-HNME  and  the  Trembly  and 
Berg  t*  determinations.  If  we  exclude  these  as  outliers,  the  difference  be- 
tween the  mean  t*  of  the  two  populations  becomes  6t*  = 0.23.  Making  a 
rough  adjustment  for  the  differences  in  the  average  epicentral  distances  re- 
presented in  the  two  groups  (about  10°)  by  multiplying  6t*  with  the  ratio  of 
travel  times  at  the  mean  epicentral  distances  the  adjusted  6t*  becomes  0.21. 
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Calculations  using  a granite  reduced  displacement  potential  for  100  kt  in 
Blandford  (1976)  shows  that  as  t*  varies  between  0.2  and  0.4,  and  assuming  no 
surface  reflection,  log(A/T)  varies  0.30  magnitude  units.  This  corresponds 
well  to  the  magnitude  differences  tabluated  by  Booth  et  al.  (1975). 

In  comparing  t*  values  for  paths  to  stations  KNOT,  DRCO  and  NLAZ  in  the 
Colorado  Plateau  province,  there  seems  to  be  an  indication  of  a directional 


EVENTS  A > 35‘ 


SHIELD  PATHS  MEAN=.22 


SHIELD  TO  BASIN  AND  RANGE  PATHS 
MEAN  = .53 


I'igure  4.  Summary  of  t*  values  for  A > 35°  , listed  in  Tables  II  and  III. 


TABLE  II 


SHIELD  AND 

STABLE  PLATFORM  PATHS 

FOR  A > 35“ 

A“ 

t* 

Central  Asia  - NORSAR 

(Noponen,  1975) 

36  - 37 

.26 

E-Kazakh  - 

- NORSAR 

(Noponen,  1975) 

38.0 

.20 

Salmon  - NPNT 

SDAC 

47.3 

.23 

N.  Zemlya 

- RKON 

SDAC 

54.0 

.15-^ 

N.  Zemlya 

- SV2QB 

SDAC 

45.8 

.18'^ 

N.  Zemlya 

- VOIO 

SDAC 

62.0 

.24-^ 

N . Zemlya 

- WFMN 

SDAC 

60.7 

.32"^ 

N . Zemlya 

- RKON 

SDAC 

54.0 

.24"^ 

N.  Zemlya 

- SV3QB 

SDAC 

45.8 

.23"^ 

N.  Zemlya 

- WNSD 

SDAC 

62.2 

.19-^ 

Average 

.22 

+ 


t 


adjusted  relative  measurements 


TABLE  III 


MIXED  (SHIELD  AND  STABLE  PLATFORM  TO  BASIN  & RANGE)  PATHS 


A® 

t* 

NTS  - NORSAR 

Frazier-Filson. 1972 

^-70 

.41 

NTS  - NORSAR 

Noponen, 1975 

-^70 

oo 

• 

Mast  - HNME 

SDAC 

''^36.6 

.95 

Piledriver  - NPNT 

SDAC 

39.9 

.36 

Knickerbocker  - SV3QB 

SDAC 

37.8 

.456 

Knickerbocker  - NPNT 

SDAC 

39.1 

.373 

Knickerbocker  - HNME 

SDAC 

36.7 

.344 

NTS  - NPNT 

Trembly-Berg 

39.0 

1.00 

N.  Zemlya  - EKNV 

SDAC 

67.4 

.57-^ 

N,  Zemlya  - GEAZ 

SDAC 

72.5 

.59"^ 

N.  Zemlya  - LCNM 

SDAC 

73.5 

.24^ 

N.  Zemlya  - LCNM 

SDAC 

73.5 

.64*^ 

Average 

.53 

+ 


adjusted  relative  measurements 
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dependence  in  t*.  P-waves  from  the  SAlilON  explosion  seem  to  be  attenuated 
while  those  coming  from  Novaya  Zemlya  are  not.  A possible  explanation  Is 
that  waves  coming  from  the  direction  of  Novaya  Zemlya  emerge  under  the  Plat- 
eau while  those  coming  from  the  Southeast  cross  part  of  the  partially  molten 
upper  mantle  region  which  is  characteristic  of  the  Basin  and  Range  and  Front 
Ranges  and  which  possibly  extends  under  the  southern  edge  of  the  plateau. 

The  Colorado  Plateau  Itself  is  characterized  by  low  heat  flow,  low  conductiv- 
ity (see  the  map  of  Gough,  Figure  20)  and,  quite  likely,  colder  and  less  at- 
tenuating mantle  (Roy  et  al.,  1972).  Table  IV  gives  a summary  of  t*  values 
for  SALMON  and  Novaya  Zemlya  shots  recorded  in  the  Colorado  Plateau.  We  have 
only  a few  values  and  the  matter  might  well  be  Investigated  further.  Spectra 
of  Novaya  Zemlya  shots  from  signals  recorded  at  stations  in  the  Basin  and 
Range  and  elsewhere  in  the  WUS  seem  to  indicate  higher  attenuation  than  for 
the  Colorado  Plateau  stations. 

To  analyze  the  Novaya  Zemlya  explosions  further,  we  have  assumed  that 
the  magnitude  vs.  t*  dependence  is  linear  and  has  the  same  slope  for  all  events 
resulting  in  an  equation  of  the  form: 

"^m  = a + bt* 
n n m 

at  station  m and  for  event  n.  We  estimated  all  a^  by  regression  analysis 
and  subtracted  them  from  the  in^'s  associated  with  each  event  to  make  the  plot 
found  in  Figure  5.  This  figure  is  similar  to  Figure  3.  The  common  slope  (b) 
was  also  estimated  and  found  to  be  -1.93  + 1.14  (95%  confidence  limits).  The 
probability  of  zero  slope,  which  would  be  equivalent  to  no  functional  depend- 
ence of  ni^*s  on  t*,  is  less  than  1%. 

The  95%  confidence  range  for  the  slope  includes  the  value  -1.50  which 
would  give  0.3  nij^  difference  for  a At*  of  0.2.  As  we  saw  above,  this  is  the 
difference  which  would  be  predicted  as  t*  varies  from  0.2  to  0.4  for  a 100 
kt  shot  in  granite  as  recorded  at  an  LRSM  station. 

Similar  differences  in  spectral  content  of  P-waves  exist  in  other  areas 


Roy,  R.  F.,  D.  D.  Blackwell,  and  E.  R.  Decker,  1972,  Continental  heat  flow. 
The  Nature  of  the  Solid  Earth  (E.  C.  Robertson,  Editor),  McGraw  Hill 
International  Series  in  the  Earth  and  Planetary  Sciences. 
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TABLE  IV 


COLORADO  PLATEAU  t*  VALUES 


From 

SALMON 

From  N. 

Zemlya 

Station 

t* 

Station 

t* 

UBO 

.592 

KNUT 

.36 

KNUT 

.639 

KNUT 

.18 

DROC 

.769 

NLAZ 

.25 

DRCO 

.29 

Average 

.66 

Average 

.27 

RELATIVE  mi, 


of  the  world.  Figure  6a  shows  a figure  from  a report  by  Noponen  showing  the 
LTM  discriminant  (normalized  third  moment  of  the  amplitude  spectrum)  at  NOR- 
SAR  for  explosions  In  various  regions  of  the  world  (black  dots) . Paths  to 
NTS  give  low  LTM  values  Indicating  the  lack  of  high  frequencies.  Note  that 
data  from  earthquakes  (open  symbols)  behave  In  the  same  way. 

Figure  6b  shows  a summary  of  t*  data  from  Table  1 as  a function  of  epl- 
central  distance.  We  excluded  the  Colorado  Plateau  measurements  from  this 
figure.  We  would  like  to  Investigate  the  possible  directional  asymmetry  In 
attenuation  In  this  region  separately  In  the  future.  The  figure  shows  a con- 
siderable scatter  In  the  t*  values.  Open  symbols  are  associated  with  paths 
which  have  one  leg  In  the  WUS.  The  rest  of  the  symbols  are  associated  with 
paths  entirely  under  the  BUS  or,  at  greater  distances,  paths  where  both  the 
ascending  and  descending  portions  are  under  stable,  shield  type  areas.  The 
distribution  of  data  with  distance  Is  primarily  determined  by  geography  of 
stations  and  nuclear  test  sites. 

Disregarding  some  extremely  high  t*  values,  and  allowing  for  some  scat- 
ter, the  existing  data  seems  to  be  In  agreement  with  the  following  general 
picture : 

At  short  eplcentral  distances  (A  < 10°)  the  two  kinds  of  paths  have  simi- 
lar t*  values,  since  all  paths  are  In  the  high  Q lithosphere.  In  the  range 
10°  < A < 20°  the  t*  values  diverge,  the  raypaths  under  the  WUS  penetrate  the 
low  Q,  low  velocity  layer  and  they  become  severely  attenuated.  At  greater 
eplcentral  distances  the  t*  values  most  likely  remain  constant,  the  difference 
between  the  two  populations  In  t*  Is  roughly  constant,  and  Is  caused  by  the 
difference  In  attenuation  along  the  WUS  leg  of  the  wavepaths.  The  minimum  t* 

difference  At*  . = .18  needed  to  explain  a regional  body  wave  magnitude  var- 

mln 

latlon  of  .25  Is  also  Indicated  In  the  figure.  Although  some  points  In  the 
two  populations  are  closer  than  that,  the  differences  In  the  mean  t*  values 
over  greater  distance  ranges,  as  shown  above,  are  sufficient  to  explain  the 
magnitude  anomalies. 

The  t*  value  .75  of  Helmberger  (1973)  used  for  computing  synthetic  sels- 
mlcs  In  the  distance  range  Indicated  Is  somewhat  higher  than  our  values,  while 


-28- 


Figure  6a.  LTM  discrimants  at  NORSAR  for  events  located  in  various  regions 
of  the  world  (after  Noponen, 1975) . Explosion  and  earthquake 
observations  are  plotted  with  filled  and  open  symbols,  respec- 
tively. (a)  Events  in  southwestern  North  America.  (b)  Explo- 
sions in  western  Russia,  earthquakes  in  southern  Greece  (circles) 
and  Turkey  (triangles).  (c)  Events  in  Central  Asia,  (d)  Earth- 
quakes at  Kuriles  (circles)  and  Sakhalin  (triangles).  Almost  all 
of  the  events  in  the  figure  are  shallow,  and  none  is  deeper  than 
100  km.  Lines  denote  LTM  theoretical  values  computed  for  t* 
values  indicated  in  the  figure  for  explosions  (solid  lines)  and 
earthquakes  (dash  lines). 


Figure  6b.  Summery  of  all  t*  data  In  Table  1 as  functions  of  eplcentral  distance 


those  derived  from  the  velocity  and  Q model  of  Archambeau^ et  al.  (1969)  seem, 
on  the  average,  to  be  on  the  low  side.  The  lowest  value  used  by  Bache,et  al. 
(1976)  t*  = .95  (they  used  values  as  high  as  1.15),  one  point  of  ours  and  the 
value  of  Trembly  and  Berg  (1968)  seem  to  be  extremely  high  relative  to  the 
rest  of  the  values. 

Losses  in  shear  deformation  would  also  imply  that  spectral  differences 
should  be  even  greater  for  short-period  S-waves.  There  are  no  t*  determina- 
tions for  short-period  S over  teleseismic  distances  in  North  America.  Sta- 
tistical analysis  of  short-period  S amplitudes  and  dominant  periods,  however, 
showed  that  the  short-period  S-waves  have  significantly  lower  amplitudes  and 
longer  dominant  periods  in  the  WUS  (Der,  Mass^,  and  Gurskl,  1975).  Figure  7 
shows  a map  of  LRSM  stations  classified  with  respect  to  S-wave  attenuation 
into  three  groups  (high,  intermediate  and  low  attenuation)  by  a statistical 
procedure.  The  region  of  high  attenuation  coincides  with  the  region  of  nega- 
tive magnitude  residuals  of  Booth,  Marshall  and  Young  (1974).  Since  the 
dominant  frequency  of  S-waves  is  lower  (.5  Hz)  in  the  EUS,  the  amplitude  ef- 
fect on  short-period  S-waves  is  not  as  great.  Figure  8a  shows  the  distribu- 
tion of  periods  of  short  period  S waves  in  the  two  parts  of  the  country.  The 
data  from  the  deep  focus  events  in  Der,  Mass$,  and  Gurskl  (1975)  were  used  to 
construct  this  figure.  The  mean  period  of  S-wave  for  each  event 

where  T_,-  and  T,_,„  are  the  averages  of  dominant  periods  of  short-period  S-waves 
EUS  WUS 

in  the  EUS  and  the  WUS  respectively,  was  subtracted  from  each  measured  value 
prior  to  plotting  a histogram.  The  histogram  show  that  the  S-wave  periods  are 
on  the  average  1 second  longer  in  the  WUS,  indicating  a significant  spectral 
difference.  The  regional  correlation  between  short-period  P and  S-wave  attenu- 
ation in  consistent  with  the  hypothesis  that  the  cause  of  regional  bias  is 
anelastic  attenuation  and  that  a considerable  part  of  the  losses  occurs  in 
shear  deformation. 


Bache,  T.  C.,  T.  G.  Barker,  N.  Rimer  and  J.  M.  Savino,  1976,  Comparison  of 
theoretical  and  observed  body  and  surface  waves  for  Kasseri,  an  explo- 
sion at  NTS,  Systems,  Science  and  Software,  SSS-R-76-2937. 

Trembly,  L.  A.,  and  J.  W.  Berg,  1968,  Seismic  source  characteristics  from 

explosion-generated  P-waves,  Bull.  Seism.  Soc,  ^. , v.  58,  p.  1833-1848. 
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REDUCED  S WAVE  PERIDD 


Figure  8a.  The  distribution  of  S-wave  periods  in  the  WS  and  EUS.  (Average 
j of  periods  for  each  event  is  subtracted). 
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From  the  above  discussion  an  internally  consistent  picture  emerges. 

The  t*  values,  determined  from  relative  spectral  ratio  measurements  correct- 
ly predict  the  absolute  magnitude  residuals  observed,  and  the  data  indicate 
that  t*  is  reciprocal  with  respect  to  source  and  receiver  positions  since 
values  of  t*  are  comparable  when  the  waves  travel  along  approximately  recip- 
rocal paths.  We  have  relatively  fewer  t*  values  for  teleselsmic  distance 
(A  > 35°),  but  the  t*  values  at  regional  distances  also  imply  about  the  same 
size  of  magnitude  differential  which  is  reasonable  since  the  teleselsmic  P- 
waves  spend  roughly  the  same  amount  of  time  in  the  low  velocity  zone  (LVZ)  as 
those  at  regional  distances.  A typical  t*  difference  of  0.2  corresponds  to  a 
magnitude  difference  of  0.3  units,  which  is  in  good  agreement  with  the  m^  dif- 
ference found  by  Booth  et  al.  (1975).  Short-period  P and  S-wave  attenuation 
have  the  same  geographical  distribution  throughout  the  United  States. 

All  of  our  available  data  presented  here  can  be  reconciled  with  two 
models;  for  shields  must  be  high  (Q^  5 1500)  throughout,  while  for  the 
WUS  the  model  is  similar  to  that  of  Archambeau  et  al.  (1969)  except  that 
in  the  low  velocity  zone  the  Is  about  one  half  of  the  value  given  by 
Archambeau  et  al. 

Figure  8b  shows  the  two  upper  mantle  Q models.  The  Q model  for  the 
shield/stable  platform  regions  of  the  world  has  a constant  Q of  1500  to  a 
depth  of  400  km.  Below  this  depth  both  models  coincide  and  Q increases  with 
depth.  This  Increase  is  required  by  the  approximate  constancy  of  t*  beyond 
A > 20®,  which  Indicates  that  Q increases  with  distance  at  the  same  rate 
as  does  travel  time.  Fine  details  of  both  models  are  uncertain,  we  are  also 
not  certain  about  the  depth  where  the  Q's  in  the  two  models  should  converge. 
The  differences  in  crustal  Q values  are  not  significant.  The  depth  distri- 
butions of  Q given  here  are  qualitatively  similar  to  those  for  Q by  Lee  and 

p 

Solomon  (1975)  who  derived  their  model  from  surface-wave  attenuation. 


Lee,  W.  B.,  and  S.  C.  Solomon,  1975,  Inversion  schemes  for  surface  wave  at- 
tenuation and  Q in  the  crust  and  the  mantle,  Geophys . R.  Astr . Soc. , 
V.  43,  p.  47-71. 
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SHIFT  OF  Mg-m^  POPULATIONS 

t 

It  was  first  observed  by  Lieberman  and  Pomeroy  (1969)  and  Basham  (1969) 
that  the  relationship  for  explosions  in  the  WUS  is  significantly  differ- 

ent  from  those  established  for  other  regions  of  the  world.  Figure  9 from 
Basham  shows  a shift  of  M^-m^  lines  for  earthquakes  (line  1)  in  the  Southwest- 
ern United  States  and  explosions  (line  5)  in  the  same  region  relative  to  earth-  * 
quake  and  explosion  lines  in  Central  Asia  (lines  4 and  7 respectively) . The 
relative  shift  in  least  squares  lines  by  various  authors  ranges  from  about  .3 
to  .5  magnitude  units  (Evernden  and  Filson,  1971).  A logical  explanation  for 
this  phenomenon  is  the  reduction  of  short-period  P-waves  amplitudes  by  anelas- 
tic  attenuation.  An  alternative  explanation,  the  enhancement  of  by  tectonic 
stress  release  seemed  to  be  ruled  out  by  ToksHz  and  Kehrer  (1972);  and  by  Masse 
(1973),  who  concluded  that  the  M^  enhancement  was  on  the  order  of  0.1  magnitude 
units  or  less.  The  magnitude  of  shift  along  the  m^  axis  seems  to  be  greater 
than  can  be  explained  by  relative  teleseismic  m^  anomalies  (Booth  et  al.,  ^ 

1975)  and  our  At*  determinations. 

i 
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Liebermann,  R.  C.  and  P.  W.  Pomeroy,  1969,  Relative  excitation  of  surface 

waves  by  earthquakes  and  underground  explosions,  J.  Geophvs.  Res. . v. 

74,  p.  1575-1590. 

Basham,  P.  W. , 1969,  Canadian  magnitudes  of  earthquakes  and  nuclear  explosions 
in  Southwestern  North  America,  Geophys.  R.  Astr . Soc. , v.  17,  p.  1- 
13. 

Evernden,  J.  and  J.  Filson,  1971,  Regional  dependence  of  surface-wave  versus 
body-wave  magnitude,  J.  Geophvs . Res. . v.  76,  p.  3303-3308. 
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SUPPORTING  EVIDENCE  FOR  HIGH  ATTENUATION  UNDER  THE  WUS 


General 

In  this  section  we  present  evidence  supporting  the  notion  of  a highly 
attenuating,  parlally  molten  upper  mantle  under  the  WUS.  Some  of  the  evid- 
ence Is  directly  related  to  short-period  attenuation  (attenuation  at  long 
periods)  and  some  Is  Indirectly  related.  Special  attention  will  be  given  to 
similarities  In  the  geographical  distributions  of  various  geophysical  varia- 
bles. 

Attenuation  of  Long-Period  Body  and  Surface  Waves 

Solomon  and  ToksHz  (1970)  published  a comprehensive  study  of  long-period 
P and  S-wave  attenuation  throughout  the  United  States  using  the  spectral  rat- 
io method.  Figures  10  and  11  show  the  regional  distribution  pattern  of  rela- 
tive t*  (<5t*)  values  obtained  for  long  period  P and  S-waves  (their  definition 
of  t*  Includes  the  multiplicative  factor  tt)  . 

It  Is  obvious  that  this  geographical  t*  distribution  Is  significantly 
different  from  that  obtained  for  short-period  waves.  Their  map  shows  a low 
attenuation  region  parallel  to  the  Pacific  Coast  and  a high  attenuation  re- 
gion In  the  northeastern  U.S.  The  regional  variation  In  t*  Is  also  greater 
than  for  short-p'^rlod  waves.  If  we  accepted  these  values  as  6t*  for  short- 
period  waves  we  would  have  considerably  more  variation  In  amplitudes  than 
that  observed.  An  obvious  way  to  resolve  the  discrepancies  Is  to  propose  the 
frequency  dependence  of  Q,  l.e.  Q would  be  required  to  Increase  with  frequen- 
cy at  different  rates  In  different  regions.  Figure  10  and  11  are  based  on 
data  from  two  deep  Peru-Brazll  border  events  which  also  may  have  led  to  dis- 
crepancies due  to  radiation  pattern  effects.  It  Is  Important,  however,  that 
the  area  of  prime  Interest  In  our  study,  the  Basin  and  Range  province  Is  char- 
acterized by  high  attenuation  for  both  long-period  P and  S-waves.  The  values 
of  6t*  for  P and  S seem  to  be  related  by  the  4:1  relationship  proper  for  los- 
ses entirely  In  shear  (Figure  12). 

Theoretical  studies  of  a solid  matrix  with  viscous  liquid  Inclusions  of 
large  aspect  ratios  show  that  such  media  attenuate  the  seismic  waves  passing 
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Figure  12.  Comparison  of  S and  P differential  attenuation  at  U.  S.  stations 
Vertical  and  horizontal  error  bars  delimit  the  range  of  calcula- 
ted values  for  each  station.  Correlation  coefficient  Is  0.67 
(after  Solomon  and  Toksoz  1970). 


through  them  (Walsh  1968,  1969).  The  attenuation  is,  in  general,  frequency 
dependent:  the  form  of  this  dependence  and  the  amount  of  attenuation  is  de- 
termined by  the  aspect  ratios,  the  properties  of  the  solid  matrix,  the  con- 
centration and  viscosity  of  liquid  inclusions.  This  model  also  requires  that 
all  of  the  attenuation  occurs  in  shear  deformation,  and  none  in  compression. 
This  seems  to  be  a realistic  model  for  a partially  molten  rock  where  some 
mineral  components  are  solid  while  others  have  melted,  and  the  shape  of  in- 
clutions,  naturally,  cannot  be  expected  to  be  equidimensional.  A solid  mat- 
rix with  inclusions  of  various  concentrations  and  shapes  can  explain  a wide 
range  of  frequency-attenuation  relationships.  The  model  implies  that  neither 
the  regional  distribution  nor  the  degree  of  attenuation  (Q)  must  be  related 
at  different  frequencies  if  the  model  changes  laterally.  Such  laterally  chang- 
ing models  can  explain  the  differences  in  the  regional  distribution  patterns 
obtained  at  different  frequencies.  It  is  significant,  however,  that  the  at- 
tenuation is  higher  in  the  WUS  for  both  short  and  long-period  waves.  The 
solid  matrix  with  liquid  inclusions  also  has  wave  velocities  lower  than  that 
of  the  homogeneous  solid  with  the  properties  of  the  matrix.  Besides,  disper- 
sion due  to  attenuation  may  make  the  delays  dependent  on  the  frequency  of 
waves.  This  model,  therefore,  also  explains  the  travel  time  delays  without 
assuming  compositional  changes.  The  sole  requirement  is  that  the  temperatures 
become  high  enough  to  melt  some  mineral  components  of  the  rock. 

Solomon  (1972)  attempted  to  model  the  observed  attenuation  and  travel 
time  delays  by  the  superposition  of  two  attenuation  mechanisms.  Although  his 
model  is  plausible,  it  is  not  unique  since  too  many  parameters  are  not  suific- 
iently  known.  Nevertheless,  detailed  studies  of  the  frequency  and  regional 
dependence  of  attenuation  and  travel  time,  as  well  as  other  data,  may  event- 
ually make  it  possible  to  build  a more  unique  geophysical  model  of  the  West- 
ern United  States. 


Walsh,  J.  B.,  1968,  Attenuation  in  partially  melted  material,  J.  Geophys. 
Res. . V.  73.  p.  2209-2216. 

Walsh,  J.  B.,  1969,  New  analysis  of  attenuation  in  partially  melted  rock, 
J..  Geophvs.  Res. . v.  74,  p.  4333-4337. 

Solomon,  D.  C. , 1972,  Seismic-wave  attenuation  and  partial  melting  in  the 
upper  mantle  of  North  America,  J.  Geophvs . Res. . v.  77,  p.  1483-1502. 
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Studies  of  long-period  surface  wave  attenuation  across  the  WUS  and  the 
EUS  (Lee  and  Solomon  1975)  also  Indicate  a low  Q layer  coincident  with  the 
low  velocity  layer  under  the  Western  United  States,  while  such  a layer  is  ab- 
sent in  the  EUS.  Figure  13  shows  their  values  of  100/Q  plotted  for  various 

p 

regions  of  the  United  States,  indicating  such  a layer  (high  100/Q  values) 

p 

under  the  WUS. 

Travel  Time  Delays 

It  has  also  been  observed  that  the  travel  times  of  body  waves  to  the  WUS 
are  consistently  late  (Cleary  and  Hales,  1966;  Doyle  and  Hales,  1967;  Hale 
and  Roberts,  1970).  Travel  times  are  not  related  directly  to  attenuation 
since  compositional  variations  can  also  cause  changes  in  travel  times  although, 
as  stated  above,  the  physical  model  of  partial  melting  of  identical  rocks  can 
explain  such  changes  without  the  need  for  assuming  compositional  changes.  Fig- 
ure 14  shows  P-wave  travel  time  residuals  from  deep  events  obtained  by  Julian 
and  Sengupta  (1973).  These  correlate  well  with  the  geographical  distribution 
of  magnitude  residuals  according  to  Booth  et  al.  and  the  attenuation  maps  of 
Der  et  al.  The  same  correlation  is  evident  on  the  map  of  S-wave  residuals  by 
Doyle  and  Hales  (1967)  (Figure  15)  and  Hales  and  Roberts  (1970)  (Figure  16); 
the  values  shown  were  derived  from  both  deep  and  shallow  events  and  thus  may 
be  contaminated  with  near  source  travel  time  effects.  The  picture  is  fairly 
consistent,  however,  showing  late  arrivals  in  the  WUS  and  early  arrivals  in 
the  EUS.  Except  for  some  mixed  delay  values  in  California  there  is  no  cor- 
relation with  the  coastal  strip  of  low  long-period  attenuation  reported  by 


Cleary,  J.,  and  A.  L.  Hales,  1966,  An  analysis  of  the  travel  times  of  P-waves 
to  North  American  stations,  in  the  distance  range  32°  to  100°,  Bull. 

Seism.  Soc.  Am. . v.  56,  p.  467-489. 

Doyle,  H.  A.,  and  A.  L.  Hales,  1967,  An  analysis  of  the  travel  times  of  S-waves 
to  North  American  Stations  in  the  distance  range  28°  to  82°,  Bull.  Seism. 
Soc.  Am. . V.  57,  p.  761-771 

Hales,  A.  L. , and  J.  Roberts,  1970,  The  travel  times  of  S and  SKS,  Bull.  Seism. 
Soc.  Am. . V.  60,  p.  461-489. 

Julian,  B.  R. , and  M.  Sengupta,  1973,  Travel  time  anomalies  for  the  global 
station  network.  Semiannual  Technical  Summary,  Seismic  Discrimination, 
ESD-TR-73-348,  Lincoln  Laboratories,  MIT,  Lexington,  Massachusetts. 
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Limits  of  100/Qg  regions  of  the  United  States  (after  Lee  and 
Solomon  1975). 


Solomon  and  ToksHz  (1970)  nor  with  their  high  attenuation  region  in  the  North- 
eastern United  States.  Comparison  of  the  magnitude  of  P and  S-wave  travel 
time  delays  shows  that  the  S-wave  delays  are  about  four  times  greater  than  the 
P-wave  delays.  This  is  consistent  with  the  notion  of  complex  shear  modulus 
caused  by  partial  melting  under  the  Western  United  States. 

Velocity  Structure 

The  crust-upper  mantle  structure  under  the  United  States  is  widely  dis- 
cussed in  the  literature.  It  is  beyond  the  scope  of  this  report  to  discuss 
all  the  relevant  information  in  detail.  Only  a few  selected  references  will 
be  quoted  in  this  summary.  Since  anelastlc  attenuation  is  generally  assumed 
to  be  associated  with  the  upper  mantle  low  velocity  zone,  only  variations  in 
the  structure  of  the  upper  mantle  will  be  discussed.  Both  body  wave  and  sur- 
face wave  studies  have  shown  that  a pronounced  low  velocity  zone  (LVZ)  is 
present  for  both  P and  S-waves  under  the  WUS.  Such  a layer  is  either  absent 
or  the  velocity  contrast  between  the  LVZ  and  the  adjoining  layers  is  much 
less  in  the  EUS.  A few  inversion  results  for  P-wave  velocity  depth  function 
are  shown  in  Figures  17  and  18  for  the  Western  and  Eastern  United  States  re- 
spectively. Figure  18  shows  a set  of  inversion  results  for  the  Western  Unit- 
ed States.  Three  models  CIT,  HOP,  HIP  and  H2P  (14a-c)  of  Archambeau  et  al. 
(1969),  the  model  SDL-UT-BRl  of  Mass4  et  al.  (1972)  (14d)  and  Johnson's  model 
CIT  204  model.  Figure  19  shows  a number  of  models  for  the  Eastern  United 
States,  a model  by  Green  and  Hales  (1968),  a superposition  of  Mass4's  (1973) 
EUS  model  (MAS)  with  a model  VIC3  by  McMechan  (1975)  based  on  the  same  data; 
and  the  model  YULKNFIO  by  Julian  (1970). 


Mass4,  R.  P.,  M.  Landlsman,  and  J.  B.  Jenkins,  (1972),  An  investigation  of  the 
upper  mantle  compresslonal  velocity  distribution  beneath  and  Basin  and 
Range  Province.  Geophys.  Astr.  Soc. , v.  30,  p.  19-36. 

Green,  R.  W.  E.,  and  A.  L.  Hales,  1968,  The  travel  times  of  P-waves  to  30®  in 
the  central  United  States  and  upper  mantle  structure.  Bull.  Seism.  Soc. 
Am. . V.  58,  p.  267-290. 

McMechan,  G.  A.,  1975,  Amplitude  constraints  and  the  inversion  of  Canadian 

Shield  Early  Rise  Explosion  data,  Bui  1.  Seism.  Soc.  Am. . v.  65,  p.  1419- 
1433. 


Julian,  B.  R. , 1970,  Regional  variations  of  upper  mantle  structure  beneath 
North  America,  Ph.D.  Thesis,  California  Institute  of  Technology. 
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Figure  18.  Selected  P-wave  velocity  models  for  the  Western  United  States 
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Figure  19.  Selected  P-wave  velocity  models  for  the  Eastern  United  States 
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We  exclude  three  models  from  our  presentation:  the  Wiggins  and  Helmber- 
ger  (1973)  models  HWA  and  HWB,  which  were  derived  by  using  data  from  a few 
close-in  stations.  The  definition  of  these  models  at  the  depth  of  the  LVZ  is 
poor  (the  authors  state  this  themselves);  the  model  HWNE  by  the  same  authors 
is  excluded  for  the  same  reason.  These  models  are  more  relevant  to  deeper 
mantle  structure,  the  450  and  700  km  velocity  gradient  regions,  that  we  are 
not  particularly  concerned  with  in  this  summary. 

Inspecting  the  rest  of  the  selected  models  it  becomes  obvious  that  the 
crust-upper  mantle  structure  is  quite  different  in  the  two  major  subdivisions 
of  the  continent.  The  WUS  is  underlain  by  an  upper  mantle  low-velocity  layer. 
In  most  inversions  one  finds  a "lid"  of  higher  velocity  just  under  the  Mohoro- 
viclc  discontinuity.  In  some  places  the  lid  may  be  missing  and  the  low-velo- 
city region  reaches  the  base  of  the  crust  (see  Model  CIT  lllP)  resulting  in 
extremely  low  velocities.  In  the  EUS  the  low  velocity  layer  is  missing 
(disregarding  an  insignificant  velocity  inversion)  the  velocities  remain  above 
8 km/sec  and  increase  with  depth  everywhere  below  the  Mohorovicic  discontin- 
uity. The  P^  (sub-Moho)  velocities  are  above  8 km/sec  everywhere  in  the  EUS 
as  shown  by  a map  of  Herrin  and  Taggart  (1962)  which  is  reproduced  in  Figure 
20.  The  high  P^  velocity  regions  in  the  WUS  do  not  necessarily  mean  the  ab- 
sence of  LVZ  or  attenuation  in  those  regions,  rather,  they  reflect  the  var- 
iations in  the  "lid"  zone.  A fairly  good  correlation  between  P^  velocity 
and  station  residuals  exists,  however,  since  low  (<8  km/sec)  P^  velocities 

occur  only  in  the  WUS.  The  WUS  is  also  characterized  by  high  S attenuation 

n 

(Molnar  and  Oliver  1969). 

Many  authors  have  derived  the  velocity  structure  for  S-waves.  The  most 
recent  comprehensive  study  of  the  subject  was  done  by  Biswas  and  Knopoff 
(1974).  They  found  a sizeable  S-wave  low-velocity  layer  under  the  WUS.  The 
same  layer  is  either  absent  or  much  reduced  in  size  in  the  EUS.  The  exis- 

Herrin,  E. , and  J.  Taggart,  1962,  Regional  variations  in  P^  velocity  and 

their  effect  on  the  location  of  epicenters.  Bull.  Seism.  Soc.  Am. , v. 

52,  p.  1037-1046. 

Biswas,  N.  N.,  and  L.  Knopoff,  1974,  The  structure  of  the  upper  mantle  under 
the  United  States  for  dispersion  of  Rayleigh  waves,  Geophys.  R. 

Astr.  Soc. , V.  36,  p.  515-540. 
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tence  of  P and  S-wave  low-velocity  layers  is  also  consistent  with  the  partial 
melting  model  described  above,  since  melting  would  reduce  both  the  wave  velo- 
cities in  the  material,  but  the  reduction  of  S-wave  velocity  would  be  greater. 
This  is  in  agreement  with  the  observation  that  the  low-velocity  layer  is  more 
pronounced  for  S-waves,  but  if  a prominent  S-wave  low-velocity  layer  is  pre- 
sent there  is  one  for  P-waves  also. 

Upper  Mantle  Resistivity  and  Heat  Flow 

Analysis  of  magnetic  variations  in  the  WUS  has  revealed  wide-spread  con- 
ductive regions  in  the  upper  mantle.  Figure  21,  a map  from  Gough  (1973)  shows 
the  areal  distribution  of  the  conductive  regions  (the  density  of  stippling  is 
proportional  to  conductivity).  The  correlation  with  short-period  amplitude 
anomalies  is  only  fair;  here  we  have  a coastal  strip  of  low  conductivity  and 
the  Wyoming  Basln-Colorado  Plateau  stand  out  as  regions  of  low  conductivity. 

The  Basin  and  Range  province  shows  high  conductivity.  The  highest  conducti- 
vities are  found  in  the  Rocky  Mountain  frontal  regions  and  a long  strip  in 
Western  Utah.  Gough  estimates  that  the  conducting  layer  is  between  50  and  100 
km  depth  in  the  upper  mantle,  but  interpretation  of  such  anomalies  is  not  uni- 
que. Vertical  topography  of  the  top  of  the  conductive  region  also  influences 
the  values  of  the  conductivity  measured.  In  areas  of  low  apparent  conductiv- 
ity the  conductive  layer  may  still  be  present  at  a greater  depth.  Correlation 
of  these  anomalous  conductivities  with  heat  flow  data  is  fairly  good  (Sass  et 
al.  1971),  as  shown  in  Figure  22.  The  heat  flow  measurements  can  be  influenc- 
ed by  many  factors,  such  as  groundwater  circulation,  sedimentary  thicknesses, 
etc. , and  as  a result  may  not  be  consistent  indicators  of  conditions  in  the  up- 
per mantle. 

Laboratory  studies  of  rocks  indicate  that  the  resistivity  of  rocks  de- 
creases by  several  orders  of  magnitude  as  the  rock  starts  to  melt.  An  in- 
crease in  temperature  of  course  will  also  lead  to  increased  heat  flow. 


Gough,  D.  I.,  1973,  The  geophysical  significance  of  geomagnetic  variation 
anomalies,  Phys.  Earth.  Planet.  Int. , v.  7,  p.  379-388. 
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Figure  22.  Heat  flow  measurements  in  the  Western  United  States 


GENERAL  COMMENTS 


Most  of  the  discussion  in  this  report  has  centered  on  the  relative  t* 
differences  along  paths  in  the  United  States.  Major  differences  in  Q struc- 

I 

ture  have  been  found,  however,  in  many  other  areas  of  the  world.  Extremely  ^ 

low-Q  regions  in  the  upper  mantle  have  been  found  behind  some  Island  arcs 
(Barazangi  et  al. , 1975;  Oliver  and  Isacks,  1967),  and  under  mid-oceanic 
ridges  (Solomon,  1973).  High-Q  characterizes  the  upper  mantle  under  conti- 
nental shields.  The  upper  mantle  under  tectonic  continental  areas  and  average 
oceans  is  probably  characterized  by  Q values  which  are  intermediate  between 
these  two  extremes.  The  major  differences  in  Q structure  found  thus  far  are 
located  in  the  upper  mantle.  It  would  be  useful  for  correct  yleld-m^  estima- 
’ tion  and  for  the  application  of  spectral  discriminants  to  compile  a worldwide 

map  of  average  upper  mantle  Q (or  t*)  values.  Knowing  the  average  attenuation 

I 

properties  of  the  upper  mantle  would  enable  us  to  correct  for  attenuation 
under  both  the  source  and  receiver.  The  same  can  be  done  to  teleselsmlc 
. body-wave  spectra.  The  deeper  mantle  seems  to  be  more  homogeneous  laterally  k 

^ in  Q,  and  it  seems  that  two  correction  factors  are  needed;  one  for  attenuation 

in  the  upper  mantle  under  the  source  for  the  downgoing  wave  and  one  for  the  * 

upgolng  wave  under  the  receiving  station.  This  suggestion  was  also  made  r 

earlier  by  P.  D.  Marshall  (personal  communication). 

Magnitude  bias  in  some  regions  can  also  be  attributed  to  causes  other 
than  attenuation.  Presence  of  thick  sedimentary  layers  under  some  regions 
(stations  along  the  Gulf  Coast  and  Florida)  is  the  probable  explanation  of 
the  high  amplitude  of  P-waves  observed  there  (see  Figure  1) . Preliminary 
calculations  for  crustal  responses  at  selected  LRSM  stations  located  both  in 


Barazangi,  M. , W.  Pennington,  and  B.  Isacks,  1975,  Global  study  of  seismic 
wave  attenuation  in  the  upper  mantle  behind  Island  arcs  using  pP-waves 
_J,  Geophvs.  Res. . v.  80,  p.  1079-1092. 

Oliver,  J.,  and  B.  Isacks,  1967,  Deep  earthquake  zones,  anomalous  structures 
in  the  upper  mantle  and  the  litosphere,  J.  Geophvs . Res. . v.  72,  p. 
4259-4275. 

Solomon,  S.,  1973,  Shear  wave  attenuation  and  melting  beneath  the  mid-Atlantic 
ridge,  J.  Geophvs . Res. , v.  78,  p.  6044-6095. 
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I i 

EUS  and  WUS  seem  to  rule  out  the  crustal  responses  as  significant  factors  in  i 

I ] 

causing  the  East-West  amplitude  anomaly  although  they  help  to  reduce  the  scat-  ! { 

I I 

ter  of  theoretical  versus  observed  magnitude  residuals  where  the  main  cause  of  I j 
magnitude  difference  is  the  difference  in  regional  t*  values.  Another  possible  ! j 
cause  can  be  focusing  of  rays  due  to  local  structure.  In  addition  to  correc- 
tions for  attenuation,  additional  station  correction  factors  may  thus  still  be  : | 

necessary. 

I 
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